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Abstract 

Toy models of interacting Pomerons with triple and quaternary Pomeron vertices in zero trans- 
verse dimension are investigated. Numerical solutions for eigenvalues and eigenfunctions of the 
corresponding Hamiltonians are obtained, providing the quantum solution for the scattering am- 
plitude in both models. The equations of motion for the Lagrangians of the theories are also 
considered and the classical solutions of the equations are found. Full two-point Green functions 
("effective" Pomeron propagator) and amplitude of diffractive dissociation process are calculated 
in the framework of RFT-0 approach. The importance of the loops contribution in the amplitude 
at different values of the model parameters is discussed as well as the difference between the models 
with and without quaternary Pomeron vertex. 
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1 Introduction 



The complete understanding of the high energy scattering processes is impossible without the calcu- 
lation of the contributions of the different unitarization corrections to the amplitude. In the QCD 
Pomeron framework, [H El E]; these corrections are arose due the self Pomeron interactions via the 
triple Pomeron interactions vertices [U [5] . These self Pomeron interactions lead to the very compli- 
cated picture of the amplitude's evolution with rapidity, see for example [6} [71 [HI [9l [TO]. 

There are few main approaches which claim that they properly describe the Pomeron self inter- 
actions. The first one is presented as a chain of the evolution equations with Balitsky-Kovchegov 
(BK) equation as a main field equation of the hierarchy, The BK equation may be formulated 
in the framework of the Color Glass Condensate Approach (CGK), |1U) . and may be properly an- 
alyzed in the terms of s-channel interacting color dipoles, [12]. This equation correctly describes 
the interaction of two non identical objects, for example such interactions as DIS on nuclei, see 
[laillliailBllITllIHllinKSO]. Another approach is the QCD Reggeon Field Theory (RFT-QCD), 
which is formulated in the momentum space and based on the standard diagrammatic calculus, 
[3 13 13 El [211 [23 [23] • in this approach the BK equation describes a resummation of the "fan" 
diagrams of the type depicted in Fig. la with only Pomerons merging vertices considered. 

Usual BK equation, describing "fan" diagrams, does not include Pomeron's splitting vertex. There- 
fore, the next natural step toward the unitarization of the amplitude is a symmetrical consideration of 
the scattering process with both vertices included. In this case the splitting vertex may be accounted 
differently in two different approximations. In the first approach the semi-classical problem of the 
interest is considered and there the sum of the diagrams depicted in Fig. lb is calculated, neglecting 
Pomeron loops contribution into amplitude. The second approach is concentrated on the solution 
of the full quantum problem basing on some effective high-energy QCD inspired models, accounting 
diagrams of Fig.lc type, see for example [24j . So far several attempts of the calculations were done 
in both directions, on the basis of QCD-RFT theory, [211 [22l [23] and in the framework of CGC and 
dipole model, [25]. The NLO and one-loop contribution into the amplitude also were calculated , see 
for example [261 [27], but the full solution still seems very far from it's completeness. In this situation 
it is very natural to consider much simpler zero transverse dimensional model, in hope that some 
important properties of solutions of this model will be faced in QCD as well. 

The RFT-0 (Reggeon Field Theory in zero transverse dimensions) model is attracted much interest 
during few last years, see [31]. The approach was formulated and studied a long time ago, even before 
the QCD era, see [281 HSl [30] . The interest to the approach is due very attractive properties of the 
model. First of all, RFT-0 is solvable at different parameters of the model. This possibility to find 
full solution for the different regimes of the theory is an important feature of the model, since we hope 
that we will understand more about RFT-QCD if we will have full and semi-classical solutions for 
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RFT-0, see examples in \32 \ \33 \ Therefore, the second important fact about the RFT-0 is that 
in this theory we can find both classical and quantum solutions for the amplitude that provides us 
with the information about the relative contribution of the loops to the amplitude. The possibility to 
consider in RFT-0 different types of the Pomeron vertices is an another property of the RFT-0 which 
is very useful and which could clarify the situation in QCD. 

In the paper we consider two RFT-0 models, in the second one together with usual triple Pomeron 
vertex we also include the quaternary Pomeron vertex. The paper is organized as follows. In the 
second section we consider the RFT-0 model with only triple Pomeron vertex. In two subsections 
of the section we introduce an whole machinery of the problem for the both quantum and classical 
solutions. In the following subsections we present the results of our calculations for the amplitude at 
different parameters of the model as well as the results for the two point Green's function ("effective" 
Pomeron propagator ) of the theory. In the Sec. 3 we solve the same problems for the RFT-0 theory 
with the quaternary Pomeron vertex included. The next section. Sec. 4, is dedicated to the possible 
application of the solution to the problem of diffractive dissociation. The Sec. 6 presents a summary 
of the results of the calculations and the last section, Sec. 7, is the conclusion of our paper. 

2 Solution of RFT in zero transverse dimension with only triple 
Pomeron coupling 

In this paper we investigate the one-dimensional problem of the interacting Pomerons described by 
Lagrangian, which in the terms of Gribov fields has the following form: 

L = - il^^ ip - fixlj"^ x/j + i Xil)^ {il;^ + i) , (1) 

where /i is the intercept of the bare Gribov field, A is the triple field interactions vertex and differ- 
entiation means differentiation on rapidity (y) variable, which is only the variable of the problem. 
Introducing the Pomeron fields, q = iip^ and p = iip we rewrite the Lagrangian in the form of real 
A coupling: 

L = qp + fiqp - \q{q + p)p ■ (2) 

The field q and p may be understood now in the terms of the Pomeron creation and Pomeron anni- 
hilation operators a: 

, d 
q^a^ = q,p^a = - — . (3) 

dq 

The Hamiltonian of the problem in the terms of the operators q and p has the following form 

H = - iiqp + Xq{q + p)p , (4) 
being the second order differential operator: 

= I + (5) 
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a) b) c) 

Figure 1: Examples of diagrams of the effective theory of interacting Pomerons with triple pomeron 
vertices: a) a fan diagram; b) a tree diagram defining the classical limit; c) a diagram with quantum 
loops. 

The full, quantum solution of the theory will coincide, therefore, with the solution of the quantum 
mechanics problem with the Hamiltonian given by EqjS) 

H^ = ^, (6) 
dy 

where the function 

n=oo 

^iy,q) = J2 Ke-^-^MQ) (7) 

n=0 

is the full quantum amplitude of the process of the scattering of two particles. Here En and (pniQ) 
denote eigenvalues and eigenfunctions of the operator Eq. O and A„ are the normalized projections 
coefficients of the eigenfunctions on the value of ^ at y = . 

The classical solution of the theory, i.e. the solution with only "tree" diagrams without loops also 
may be obtained in the given framework. It is simply the solution of equations of motion for the 
Lagrangian given by Eq. [2] 

2.1 The quantum solution of the model 

An approximate solution of the Eqj6]was found many years ago, see [281 [291 [30] . Analytical solutions of 
the Eqj6| at large value of the ratio fi / X = g for ground state and first exited state were considered 
at the same time as well. Nevertheless, the spectrum of the theory and eigenfunctions for next 
exited states at large value of g as well as the solutions for arbitrary values of g were not obtained. 
Therefore, we want to define the procedure for the numerical calculations of the spectrum of the theory 
and corresponding eigenfunctions for the cases of the strong triple Pomeron coupling (small g ) and 
week coupling (large g ) as well. In the following we will be more interesting in the case of the small 
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triple coupling as the case which is mostly correspond to the real perturbative QCD situation. So, we 
need to solve the second order differential equation (see more detailed in [29]): 

dq dq^ Xq dy 

with the initial and boundary conditions on the function ^(y, 



[Q - q) — 7z — + — 1Z2 — = — — 7— ' (8) 



^{y = o,q) = E Mq) = i{q) (9) 

n=0 

^{y,q ^ 0) oc g (10) 
"^{y,q 00) oc const. , (11) 

where form of I{q) depends on the particularly considered physical problem. Now it is constructive 
to make a transformation of the eigenfunctions: 

Uq) = e^'-'^" fn{q) . (12) 

In this case the term which is proportional to the first derivative over q is eliminated leading to hermi- 
tian Hamiltonian for the problem. After this transformation we obtain a Shredinger type eigenvalue 
equation: 

+ ^-1^1-6) Uq) = -—/n{q) (13) 
with the following boundary conditions on the function fn{q)'- 

Uq ^ 0) OC g (14) 
/„(g ^ 00) oc e-^'/^ + ^^-Z^ _ (15) 

Our problem of interest is the scattering of two particles, therefore we define initial condition for the 
^{y = 0, q) in the eikonal form (see more in p9]) : 

^{y = 0, (?) = I{q) = I{q , q^^t) = 1 - e"^^-* , (16) 

here qext is a value of the source for the Pomeron field at zero rapidity. The function ^{y,q) results 
the evolution of interacting Pomeron fields for the value of the Pomeron-external particle vertex equal 
to qext at zero rapidity, till the value of the single Pomeron field equal to q at rapidity y. 

The last ingredients of the theory are projection coefficients A„ of the eigenfunctions of the solution 
on the initial state I(q). Having in mind, that our eigenfunctions are orthogonal on interval of q from 
to 00 with weight function F]y{q,Q) : 

/ fn{q) fmiq) Fw{q, g) dq = 5nm , (17) 
Jo 



we obtain for X„ : 



. I ^ _ /o°° fn{q) I{q , qext) Fwjq, q) dq 

!^ fl{q)Fw{q.Q)dq ' ^^'^ 
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where the weight function Fw{q, g) has the following form (see |29] and also [?]): 

Fw{q,Q) = . (19) 

q 

The numerical solution of Eq|T3] with the boundary conditions given by Eq|T3] is the solution for 
the eigenvalues and eigenfunctions of the second order differential equation with two boundary value 
conditions. These solutions for eigenfunctions and eigenvalues may be found for each values of fi 
and A, i.e. for the different values of parameter g . The obtained values of eigenfunctions En and 
eigenfunctions fniq) allow to complete the calculation of the full quantum scattering amplitude and 
, therefore, solve our problem. So, considering the scattering of two particles , where the vertex of 
interaction of the first particle with Pomeron is equal qi and the vertex of interaction of the second 
particle with Pomeron is equal q2, we define the full , quantum solution for the scattering amplitude 
at given rapidity y as 

n=oo 

Hy,Q = q2) = E A„(gi)e-^"S'e(''2-^)'/Vn(92) , (20) 

n=0 

with the XniQi) given by Eq. [18] for the known values of eigenfunctions, see Fig.lc. The solution for 
the amplitude, Eq. [201 does not depend from which value of Pomeron field, qi or q2, the evolution 
begins. If we consider the evolution of the Pomeron from q2 to qi we will obtain the same answer for 
the amplitude, as for the case of evolution from qi to q2- 

2.2 The classical solution of the model 

In order to see the role of the Pomeron loops in the scattering amplitude it is important to calculate 
the classical solution of our problem, given by the "net" diagrams of Fig. lb. The comparison between 
two solutions will show a relative weight of the loops in the scattering amplitude and, therefore, an 
applicability of the classical solution at given values of the theory parameters. The required classical 
solution we obtain solving equations of motion for the Lagrangian given by Eq[2| with the sources of 
Pomeron fields at zero rapidity and final rapidity of the process Y: 

L = ^qp - ^qp + f^qp - Xq{q + p)p + q{y) po{y) + qo{y) p{y) , (21) 

where poiy) qo{y) are the sources of the Pomeron fields. The equations of the motion for the fields p 
and q are the following: 

q = fiq - Xq^ - 2Xqp (22) 
p = - + Ap^ + 2Xqp (23) 
qo{y) = qi 5{y) (24) 
po{y) = q2 5{y-Y) . (25) 



6 



This system of equations is another example of two value boundary problem for the system of first 
order differential equations, and may be also solved for different values of parameters /x, A, qi, q2 ■ 
The solutions of the system Eql22] at given value of final rapidity Y and given values of parameters 
fi, A, Qi q2 we will denote by {qdy) , Pc{y)} ■ With the solution {qdu) , Pc{y)} the "net", classical 
amplitude , represented by diagrams Fig. lb, is defined by standard way as: 

^ciY,qc,Pc) = 1 - e-^^^'"-^^) , (26) 

where 

SiY,q„pc)= r L{Y,qdy).Pc{y))dy = ]-{qipM+qc{Y)q2) + ^ C {ql{y)Pc{y)+qc{y)pl{y))dy . (27) 

The amplitude EqJ2S] describes the eikonalized interactions of "net" diagrams and for symmetric 
boundary conditions, qi = q2, the amplitude is invariant under the duality transformations: 

p q and y ^ Y - y . (28) 

The interesting feature of the classical solution of the system defined by the Lagrangian Eql^ 
is that starting from some critical rapidity Yc the solution of equation of motion is not unique for 
Qi = Q2 < Q- From the rapidity Y^, there are at least three classical trajectories {^c(y)iPc(y)} which 
locally minimize the action and with increasing of the rapidities the number of the trajectories is 
growing. Each from them provides the local minimum of the action and the amplitude of the theory, 
therefore, must be rewritten in the following form: 

^dY) = 1 - 5: A, exp{-5(y, q'MJM)}^ (29) 

i 

where Aj is the quantum weight of the corresponding classical trajectories. In the following consid- 
eration these weights we take equal 1 or -1, depending on the type of trajectory, see more detailed 
consideration of this problem in [3^ . 

In our calculations of the classical amplitude for the symmetrical case qi = q2, we always will use 
three solutions of the equations of motion. The first one, {ql{y),pl{y)}, is the symmetrical in the 
sence of Eql28| which trajectory will be similar to the trajectory number 1 in the Figl2j Another 
two dominant solutions, which we use in the calculations and which arise from some rapidity Yc , are 
the solutions for the trajectories similar to the trajectories 2 and 3 in the Figl2i Separately each of 
them is not symmetrical in the sense of the symmetry transformation given by the Eg 1281 but instead, 
under the duality transformation of rapidity given by the Eq j28| there is a pair symmetry of these two 
solutions 

Qciy) = Pl{Y-y) and pl{y) = ql{Y-y). (30) 

It is interesting to note, that asymptotic behavior of these two classical solutions is well known, each 
of them may be described by "fan" diagram amplitude of Fig. la, see |29j, and being the dominant 
contributions these solutions lead to the "fan" dominance effect, see |32j . 
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Figure 2: Classical solutions of the RFT-0 with only triple Pomeron vertex: the {q.,p\ trajectories 
obtained for Y = b > Yc , q = p = 0.7, p = 5. 

The final solution of our problem in the classical approximation we can write with the help of these 
three particular solutions as 



^,{Y) = l + exp{-S{Y, ql{y),pl{y))}-exp{-S{Y, ql{y) ,pl{ym-exp{-S{Y, ql{y) ,pl{y))} , (31) 



which is symmetrical in respect to the duality transformation EqJ28l In our plots for the symmetrical 
interactions, qi = q2, the classical solution will be always presented by amplitude of Eq l31l which 
contains three parts beginning from critical rapidity Yc, and only one part, symmetrical solution, at 
rapidities smaller than Yc- In the asymmetrical case the symmetry is broken initially, and, therefore, 
only one "fan" configuration survives in the classical solution. 

It is important to underline, that this picture for classical solution holds for each from the consid- 
ered models, for RFT-0 with triple Pomeron vertex only and for RFT-0 model with both triple and 
quaternary vertexes. 

2.3 Parameters of the model 

First of all, we define the range and value of parameters for which the calculations were done in 
both cases, quantum and classical. We calculated quantum and classical amplitudes, ^{Y, q) and 
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Table 1: The eigenvalues of Eq[T3l for different values of parameter g. 



^c(^> Qc, Pc) 1 for the following value of g = fi / X: 

g = 5 (// = 0.2, A = 0.04). (32) 

The reason of the main attention on this value of g is very simple. The coupling constant in RFT-0 
may be defined as = 1 / g and at = 5 we have, therefore, = 0.2 . This value for the 
coupling constant is reasonable also in QCD that provides the possibility of some analogies between 
two theories. The value of the external sources we take as follow 

• symmetrical case: 

gi = <72 = 0.1 - 1 ; (33) 

• non-symmetrical case: 

qi = 0.2 ; ^2 = 0.3 - 1 . (34) 

As an example of the case of the strong coupling , i.e. large value of the triple vertex, we also will 
present the quantum solution for the symmetrical case of interaction at 

^ = 1 (^ = 0.1, A = 0.1). (35) 

Now it is instructive to consider the spectrum of the RFT-0 theory at different g. The Tabl^ 
presents the found eigenvalues of the theory, for g = 5 as well as for g = 1 and for g = 3 . The 
asymptotic behavior of the amplitude is clearly seen from the Figj3] which represents the behavior of 
the ground state as a function of g . Being important at small rapidity, all exited states, i.e. states 
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Figure 3: The ground state of RFT as a function of g. 

with eigenvalues Ei, i = I..00 , at larger rapidities are decreasing very rapidly to zero. Therefore, at 
asymptotically large rapidity only the ground state survives, which nonzero origin is due the tunneling 
effect of the order 1'^ which arises between the Coulomb and harmonic wells of the potential of 
the EqUni see [281 [291 130] . At large enough rapidity the amplitude for any value of q is zero: for small 
Q it happens rapidly already at small rapidities, at large q huge rapidity needs. As we will see further, 
this is a particular property of the model with only triple Pomeron vertex, the theory with additional 
quaternary vertex has precisely zero ground state. 

2.4 Numerical results: loops versus "net" diagrams 

The results of the calculations of the scattering amplitude iT for the quantum and classical cases are 
presented in the FigSj-FiglTj for the symmetrical values of the external sources. In the Figj8]-Figl9j 
the results for the non symmetrical case are represented. These results are presented in the form of 
3D plots as well, as in the form of contour plots. The white color of the contour plot corresponds 
to the maximum value of the z axis function of the 3D plot, and the black color , correspondingly, 
to the minimum value of the same function. The amplitude is given as a function of two variables, 
rapidity and value of sources. In a symmetrical case the sources are the same, in non symmetrical 
case we fix one sources, q\ = 0.2 , and the second variable is the value of the second, external source 
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q2 = 0.3 — 1 . We also present the absolute ratio of the difference between quantum amplitude and 
classical amplitude to the quantum amplitude, \^{Y, q) — ^c(X^ Qc, Pc) \ / ^(^i q), in order to illustrate 
the relative contribution of the loops in the full solution comparatively to the classical solution. The 
rapidity variable Y of the plots is scaled, i.e. in the plots the variable fiy was used instead usual 
rapidity y. Therefore, in order to obtain the value of the amplitude at some physical rapidity it needs 
to divide the Y variable of the plot on intercept fi : Yp^yg = Ypi^t I • For any value of intercept [i 
we have the plot of the amplitude for rapidities Yphys = Ypiot / at triple Pomeron vertex A equal to 
A = n / Q . We see, that for the maximum value of the l^^ot = 8 at intercept, for example, ^ = 0.3 
the normal rapidity is large: Ypiiyg ~ 27 and covers the reasonable rapidity range which may be 
interesting in the high energy physics. 




Fig. gD-a Fig. g])-b 

Figure 4: The quantum amplitude '^{Y, q) of the 3P vertex model in the form of 3D and contour plots 
at Q = 5 as a functions of scaled rapidity Y and symmetrical values of external sources qi = q2- 

We also present a quantum amplitude calculated for the small value of ^ = 1. This value of g 
means large value of triple Pomeron vertex and relatively large value of the ground state energy, see 
tabled) From the plot FiefTOlis clearly seen, that at large enough rapidities the amplitude approaches 
zero as it must be for such value of g . 

2.5 "Effective" Pomeron propagator 

With the knowledge of the spectrum of RFT, the calculation of the "effective" Pomeron propagator 
P(.ff{y,q) is easy task. First of all, we change the initial condition for our equations, instead EqfT6l 
now we have 
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Fig. ©-a 



Fig. dSD-b 



Figure 5: The classical amplitude ^c(X, Qc, Pc) of the 3P vertex model in the form of 3D and contour 
plots at g = 5 and presented as a functions of scaled rapidity Y and symmetrical values of external 
sources qi = q2- 




Figure 6: The ratio \^{Y, q) — ^dY, qc, Pc) \ / ^{Y, q) in the 3P vertex model is represented by 3D 
plot and contour plot. The plots are presented as a functions of scaled rapidity Y at g = 5 and at 
qi=q2 = 0.2 - 1. 
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Fig. dZD-a 



Fig. dZD-b 



Figure 7: The same ratio \^{Y, q) — ^ciX^ Qc, Pc) \ / ^{Y, q) in the 3P vertex model is represented by 
3D plots, now at the values of the sources qi = q2 = 0.1 — 0.4 in the Fig.7-a and at the values of the 
sources qi = q2 = 0.4 — 1 in the Fig.7-b correspondingly. The plots are presented as a functions of 
scaled rapidity Y at g = 5 . 




Fig. dSD-a Fig. dSD-b 

Figure 8: The quantum amplitude ^{Y, q) of the 3P vertex model for the non symmetrical case is 
represented by 3D plot and contour plot. The plots are presented as a functions of scaled rapidity Y 
at g = 5 and at fixed qi = 0.2 with q2 = 0.3 — 1. 
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Fig. ©-a 



Fig. ©-b 



Figure 9: The ratio \^{Y, q) — ^dY, Qc, Pc) \ / ^{Y, q) in the 3P vertex model is represented by 3D 
plot and contour plot at the nan symmetrical values of the sources: qi = 0.2 is fixed and q2 = 0.3 — 1 
is a variable. The plots are presented as a functions of scaled rapidity Y at g = 5 . 




Y Y 
Fig. (UnD-a Fig. l^ -h 

Figure 10: The quantum amplitude ^{Y, q) of the 3P vertex model in the form of 3D and contour 
plots at Q = 1 as a functions of scaled rapidity Y and symmetrical values of external sources qi = q2- 
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^(y = 0,g) = I{q) = I{q,qext) = qQext- (36) 

The scattering amplitude ^{y,q) now describes the Green's function of transition from one Pomeron 
created from source q at rapidity zero to any number N Pomerons interacting with sources qext at final 
rapidity. 
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Fig. (fTD-b 



Figure 11: The Pampiy ■, q) and Peff{y,q) functions of the 3P model at g = 5 andqi = q2 =l/g =0.2 
as a functions of scaled rapidity Y . 

Therefore, in order to obtain requested "effective" Pomeron propagator with only one Pomeron 
interacting with one source at final rapidity, we also need to take the derivative of ^{y, q) over Pomeron 
field q at q = and multiply obtained function on Pomeron field q : 

■d^{y,qy 



Peff{y,q) 



■ Q 



(37) 



dq / q=0 

The Peff{y, q) function is simply the second term of the Tailor expansion of the ^{y, q) around q = 0. 
Dividing obtained propagator Pf,ff{y^q) on the sources qqext, we obtain the Green's function of the 
theory, which does not depend on the values of the sources: 

Peff{y,q) 



Pamp{y,q) 



qqext 



(38) 



Pampiy,q) does not depend on the sources of the Pomerons fields and in order to obtain requested 
propagator with given sources we need simply to multiply Pampiy,q) on the values of these sources 
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qqext- FiefTT] presents a plot of the Pampiy,Q) function and plot for the Peff{y,q) function at q = 
qext = ^/q = 0.2. The interesting question, which we can address in this calculations, it is the 
question about the importance of Pf.ff{y-, q) in RFT-0. Namely, let's define with the help of P(,ff{y, q) 
the eikonalized "effective" Pomeron amplitude 

Kffiy^Q) = 1 - e-^^//(^''^) (39) 

which account loops in the "effective" propagator and does not account interactions between different 
propagators and between propagators and external sources. This amplitude could clarify the precision 
of this approximation in comparison with the full solution. So, the required ratio is presented in the 
Fi^. 




Fig. (IH-a Fig. ^-h 

Figure 12: The ratio q) — ^f^j-{y, q) \ / '^{Y, q) in the 3P model in the form of 3D and contour 

plotshe 3P model at g = 5 as a functions of scaled rapidity Y and symmetrical values of external 
sources qi = q2 = 0.1 — 0.2. 

3 Solution of RFT in zero transverse dimension with triple and 
quaternary Pomeron vertices 

In this section we discuss a RFT-0 model given by the following Lagrangian: 

L = qp + fiqp — Xq{q + p)p + A q^ p^ , (40) 



16 



where new quaternary Pomeron vertex A' is introduced in comparison to the EqlJl In the terms of 
the q, p operators, see EqiSj the Hamiltonian of the problem has the form 

H={^,q- W)^^ + {M - A'?')^ , (41) 

and the full, quantum solution of the theory will coincide with the solution of the quantum mechanical 
problem determined by the new Hamiltonian Eql41t 

i?^ = ^, (42) 
dy 

where the function 

n=oo 

*(y,g) = Ke-^^-^MQ) (43) 

n=0 

is the full quantum amplitude of two scattering particles. As before. En and (j)n{q) denote eigen- 
values and eigenfunctions of the Hamiltonian given by the Eq l41l coefficients A„ are the normalized 
projections coefficients of the eigenfunctions on the value of ^' at y = . 

The value of the quaternary vertex A' is important for our following calculations. If the vertex 
is very small, A < \ / Q then all results of the Section 1 will stay the same with only very small 
corrections to the amplitude, it is was known many years ago, see [29]. Our calculations support this 
conclusion. Therefore, for our calculations more interesting to take another value of the vertex, so 
called "magic" value of A' , which is given by the following expression 

A' = - . (44) 
Q 

The important property of the theory with this "magic" value of A' is that in this theory , unlike 
previous cases, the ground state is precisely zero, as we will see later. Another important feature of the 
Hamiltonian Eq2T]with the "magic" value of the vertex is that this model has precise correspondence 
with the s-channel reaction-diffusion models, see ||35j) and with conformal approach for the QCD 
Pomeron, see [5^ [57]. 

As in the previous model, in this section we will also consider only the calculations with the value 
of the parameter q = b. Nevertheless, the qualitative understanding of the behavior of the amplitude 
at smaller values of q also will be clear. 

3.1 The quantum solution of the second model 

The quantum solution of the Hamiltonian EqUT] is a solution of the following second order differential 
equation 

/ , , <f-^{y,q) , , , , d-^{y,q) d^{y,q) 

A, (A/A -^)^^+ Ag(^>-g)^^ = ^^ (45) 
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with the initial and boundary conditions on the function ^{y, q): 

n=oo 

^{y = 0,q) = Mq) = Hq) (46) 

n=0 

^{y,q ^ 0) <x q (47) 
^(2/) q ^ g) oc const. , (48) 

where as before, the form of the function I{q) depends on the particular physical problem. Very 
important property of Eqj35]is the existing of the ground state (j)o{q) with the zero energy £'0 = . 
In this case the solution for EqJ45lis trivial: 



9o 



(q) = 1 - (49) 



This ground state is not orthogonal to the other eigenfunctions of the Hamiltonian, (pi, i = I..00, 
which are orthogonal to each other: 

re 







with the weight function 



where 



4>i{q) (j)j{q) Fw{q, q) dq = Const 6ij (50) 



Fw{q,Q) = —, r- (51) 

(i{q - q) 



(Pi(q = 0) = 0, i= I..00; (52) 
"^i(9 = £•) = 0, i = 1 .. 00 . (53) 

Using these properties of th eigenfunctions we obtain for the Aq projection coefficient: 

, Hq) Hg,qext) .... 

4>o{q) 1 - e s 



All other projection coefficients, Aj, i = 1 .. cxo , may be calculated using the following expression: 

. f . _ lo 4>i{q) I{q , Qext) Fwjq, g) dq _ 4>i{q) 4>Q{q)Fw{q, q) dq 

Io'<t>K<l)Fw{q,g)dq Io'4>Uq)Fw{q,g)dq ' ^''^ 

For the functions (pi, i = 1 . . 00 the equation EqjlSl has the form 

d'^4>i{q) , d(j)i{q) Ei 

—p^ + Q—, = -Tw \(t>i{q)^ (56) 

dq'^ dq ^ q[g — q) 

and with the help of the transformation 

Uq) = h{q)e-'^'^ (57) 

this equation gets the following standard Shredinger form: 

d^ f^{q) f( ^ Ei 

-^-T^^^'^ = - X'qig-q) ^^^'^ 
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for the fi{q) functions. The function fi{q) are defined on the edges of interval as 

/,((? = 0) =0, i= I..00; (59) 

Mq = Q) = 0, i= I..00, (60) 

and they are orthogonal each to other on the interval from to g with the new weight function 

Our problem we also can solve using hermitian Hamiltonian defined by Eql58] with the projection 
coefficients 

^ lo' n{q)Iiq,qe.t) e<^<i/^Fl{q,g)dq _ ^ Jj hiq) Mq) e^<i / ^ F^^jq, g) dq 

lo' fHq)FUq^8)dq ' Si fKq)FUq,e)dq 

As was done in the previous section, solving the two value boundary problem for differential equa- 
tion EqlSU] or EqJ58| we obtain the spectrum and eigenfunctions of the model and find an amplitude 
for the different values of the sources gi, q2 and parameter g: 

^{y,q = q2) = ^^'''t^^^"/"""^ + Uqi)e-''-yK{q2) ■ (63) 

1 — e~s ^ 

n=l 

3.2 The classical solution for the second model 

Comparison between the contribution to the amplitude of diagrams with loops and without loops may 
be performed again if we additionally will calculate classical solutions for the Lagrangian Eql40l 

L = ^qp - ^qp + nqp - Xq{q + p)p + X' q^ p'^ + q{y)po{y) + qo{y)p{y) ■ (64) 
The equation of motion for the fields p and q in this case are the following: 

q = fiq - Xq^ - 2Xqp + 2X' qp^ (65) 

p = —fJ-p + Xp^ + 2Xqp — 2X q^ p (66) 

qo{y) = qiKy) (67) 
po{y) = q2Ky-Y) . (68) 

The solutions of these equations, of course, are different from the solutions of the equations of motion 
for the Lagrangian Eq l21l but the representation of the amplitude ^ciX) in the terms of three classical 
solutions of the system Eql65]is the same as in the Subsection 2.2 ( see also Fig JT3]) : 

^,{Y) = l + exp{-S{Y, ql{y),pl{y))]-exp{-S{Y, ql{y) ,pi{y))] - exp{-S{Y, ql{y) ,pl{y))] . (69) 

Therefore, here we will not repeat all steps of the Subsection 2.2. The picture for the classical solution 
derived for the model with only triple Pomeron vertex is the same for the model with both triple 
and quaternary vertices. The only difference between two models is the value of critical rapidity 
Yc from which the additional solutions arising, breaking the target-projectile symmetry of the initial 
symmetrical solution, but for our consideration it is not important. 
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Figure 13: Classical solutions of the RFT-0 with triple and quaternary Pomeron vertexes: the {q^p} 
trajectories obtained for y = 5 > Yc, q = p = 0.7, p = 5. 

3.3 Parameters and asymptotic behavior of the amphtude in the second model 

As in the previous case, for the second model we will consider quantum amplitude ^{Y, q) for the 
following value g = /U / A : 



g = 5 (/i = 0.2, A = 0.04); 
and for the following values of the external sources (only symmetrical case ): 

gi = Q'2 = 0.1 - 1. 



(70) 



(71) 



The ground state of this model has precisely zero eigenvalue, as we showed before, see Eqll9] and 
for other eigenvalues of the model see Table Indeed, now, at asymptotically large rapidity, the 
amplitude is fully defined by the ground state: 

^asymp.(q) = >^oMQ) = Z "3 • (72) 

1 — e ^ 

In the case of the eikonal type initial conditions I{q,qext) the quantum amplitude ^'(y, g) at q = g 
has the following form 

/(£),(?ext) = 1 - e-^^-S (73) 
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Eo 


El 


E2 


E3 


E4 


E5 








i5^9 


g = 5 





0.182 


0.189 


0.307 


0.338 


0.42 


0.513 


0.626 


0.756 


0.905 



Table 2: T/ie eigenvalues of EqlSEfor the value of the parameter g equal 5. 



and, therefore, we obtain the following asymptotic behavior of our quantum amplitude : 

"fasymp.iq) = ^^72 (74) 

1 — e ^ 

For the "effective" Pomeron propagator, Peff{y,q) the initial condition is different: 

I{0,Qext) = QQext- (75) 

Using the definition of the propagator given by EqJST] at asymptotically large rapidity we obtain: 

p;;j"^p-(g) = (gg^^^Hgg) . (76) 

As we see, for each particular values of sources qi,q2 and parameter g neither quantum amplitude 
nor full Pomeron propagator are not decreasing to zero at asymptotically large rapidity, approaching 
instead some constant values defined by EqlTH and EqjTB] correspondingly. 



3.4 Numerical results for the the second model 

The results of our calculations of the scattering amplitude iT for the quantum case are presented in 
the same from as for the first model. Due the similarities between these two models, in this section 
we present only quantum solution for the amplitude for the symmetrical case of interactions, see 
Figll4[ We also present the ratio between the difference of the triple Pomeron vertex model quantum 
amplitude and quantum amplitude of the second model to the the quantum amplitude of the first 
model : | ^'^^ - ^-^^ | / ^'^^ , see FigHH 

3.5 "Effective" Pomeron propagator in the second model 

The "effective" Pomeron propagator we define as in the previous model 
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Fig. (HlD-a 



Fig. (Hl-b 



Figure 14: The quantum amplitude ^{Y, q) of the second model in the form of 3D and contour plots 
at g = 5 as a functions of scaled rapidity Y and symmetrical values of external sources qi = q2- 




Figure 15: The ratio of the quantum amplitudes \ — \ / in the second model in the form 
of 3D plot at Q = 5 as a functions of scaled rapidity Y and symmetrical values of external sources 



22 



with the initial condition given by 



^(y = 0,g) = I{q) = Iiq,qext) = QQext 



(78) 



The asymptotic behavior of the P(,ff{y, q) at asymptotically large rapidities is defined by expression 
EqiTSl 



(QQext) (Qq) 



1 — e ^ 

For the values qext = q = 1 / Q , we obtain that at large rapidity this function reaches unitarity limit 
see plot FigJ16l and more detailed derivation in [35j. The plot of the Green's function of the theory 



is presented in FigfTGl 



amp 



Peff{y,q) 

q qext 



(80) 
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Fig. (HED-b 



Figure 16: The Pamp{y,q) md Peff{y,q) functions of the second model at g = 5 and qi 
1/g = 0.2 as a functions of scaled rapidity Y . 
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4 DifFractive dissociation process in RFT-0 



As an example of applications of the methods developed in the paper, we present calculations of the 
diffractive dissociation process in the RFT-0 with only triple Pomeron vertex. The algorithm of the 
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Figure 17: The single diffractive dissociation process for rapidity Yi at total rapidity Y . 
calculations is defined as follows. We have the quantum solution of the theory 

n=oo 

= A„e-^"J^0„(g) (81) 

n=0 

defined at all range of rapidity Y at some values of the sources qi and q2- Let's now consider this 
solution at the rapidity interval Yi < Y. With the help of the ^i{Yi,q) we define another function 

n=oo 

My,q) = E A„e-^"(^-^^)<A„(g) (82) 

n=0 

on the rapidity interval Yi < y < Y with the following initial condition: 

n=oo 

^2iy = Y,,q) = E A„ Mq) = (*i(Y'i,g))2 . (83) 

n=0 

The illustration of this construction see in the FigJ17[ 

The only numbers that we need to calculate now are the coefficients An • But their calculation is 
trivial, they simply 

I /o°° Mq) ( MY,,q) f e-fa-^)^ Fw{q, q) dq 

" ^ ~ lo^ M<D FwiQ, e) e-(^-^)^ /2 dq ' ^''^ 

where the weight function Fw{q, g) is the same as in the EqJ19l The coefficients A„ fully determine 
requested ^2(11, q) function, which represents the differential cross section of the single diffraction 
process at given and fixed rapidity interval 

y = Y2 = Y-Yi 
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Integrating this function over rapidity interval li < y < Y we obtain as the answer the sum of the 
total diffractive dissociation cross sections on this rapidity interval and elastic cross section 

'^Id + ^ei = r My,Q)dy (85) 
So, the value of the single difractive cross section integrated over rapidity interval Y2 is the following 

crsD= r My,q)dy -{^iiY,q)f , (86) 

where Y is the total rapidity of the process, Yi is the rapidity "gap" of the process, I2 is the value 
of the rapidity taken by produced diffractive state and elastic cross section of the process is defined 
as (Te/ = (^i(y, g))^. In general we do not expect that the numbers given by RFT-0 theory will be 
correct. The only interesting value in the given model, therefore, is the ratio of the single diffractive 
cross section to the total cross section. 

Tot 

R(Yi) = ^ (87) 

where atot = '^^iiY,q). Considering this ratio we can not prove, of course, that the same value of 
the R we will obtain in QCD as well, but, nevertheless, it is interesting to see what new about this 
ratio we obtain using full quantum solution for the amplitudes. Another interesting ratio, which we 
can calculate in the given framework, this is a ratio of differential single diffraction cross section to 
the total cros section: 

^2{Y2) 



Ri{Yi) 

O'tot 

In i?i ratio the dependence on diffractive mass Y2 is present, and it allows to as to trace the energy 
dependence of the ratio of the diffractive state to the total cross section. Therefore, the following 
calculations are presented. We fix the value of diffractive state Y2 equal 1-4 units of scaled rapidity as 
a one variable and as a second variable in the 3D plot we consider the value of total rapidity Y, which 
we take equal 4-8 units. These plots we make for the two different symmetrical values of the external 
sources: 0.2 and 0.7. Obtained plots, see FigJ18lfor the ratio R and FigJlQIfor the ratio show the 
changes of the ration accordingly to the total rapidity changes. 



5 Discussion of results 



The are the following main results of the paper. First of all, we obtained the full quantum solutions 
for the models and compared it with the classical solutions, illustrating the relative contribution of 
the loops to the amplitude. In order to find an useful approximation to the full quantum solution we 
also considered an eikonalized amplitude which was built with the help of the full two point Green's 
functions. All these amplitudes were calculated for different values of external vertexes that allows 
to trace the applicability of each approximation separately in different parameter's regions. Another 
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Fig. ([IHD-a Fig. (fT8|)-b 




Fig. (HHD-c Fig. ([HD-d 

Figure 18: The ratio R for the case of symmetrical sources qi = q2 = 0.2 in the Fig. il8\) -a - Fig. il8\}- 
b and for the symmetrical sources qi = q2 = 0.7 in the Fig. U8\) -c - Fig. Iil8\)-d correspondingly. The 
ratios are plotted at g = 5 as a functions of scaled rapidity Y and scaled rapidity variable Y2. 
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Fig. dED-a 



Fig. dED-b 



Y 4 p=5 




4 

Y_2 Y 

Fig. ([l9])-c Fig. (fT9|)-d 

Figure 19: The ratio Ri for the case of symmetrical sources qi = q2 = 0.2 in the Fig. U9\) -a 

- Fig. il9\) -b and for the symmetrical sources qi = q2 = 0.7 in the Fig. il9\)-c - Fig. U9\) -d 
correspondingly. The ratios are plotted at g = 5 as a functions of scaled rapidity Y and scaled 
rapidity variable Y2. 
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important result is a comparison of quantum solutions with and without quaternary Pomeron vertex 
that clarify the contribution of this vertex to the process amplitude . 

Plots Fig|l]-Figl2] represents the result of the calculations in the RFT-0 approach with triple 
Pomeron vertex only. The quantum solution of the model is presented in the plot FigjU whereas 
the classical solution is presented in the plot Figl5l The comparison between these two solutions could 
be found in the plots Figi6]-Figi7J From the FigU] we see, that our full solution for small values of 
external sources does not reach unitarity really. Instead of "black" disk the "grey" one is achieved. 
This is may be easily explained if we will remember that our initial conditions have a eikonal form, 
whereas in the s-channel model, see [35j and references therein, the propagator form of initial condi- 
tions is always applied. In RFT-0 model the black disk limit, as it must be, achieved only at large 
values of the sources, that corresponds to the nuclei interaction in usual QCD, see [211 [22l [38] . Other 
interesting results are represented in the Fig[6]-Fig[7| Indeed, the natural question, which in fact is 
very important in any practical calculations with BK equation involved, is a question about applica- 
bility of classical solutions of the model. Namely, at which values of external sources the difference 
between the classical and quantum solutions is acceptable from the point of view of the precision of 
the amplitude calculations. In real QCD this answer may be obtained only approximately, at least. In 
RFT-0 calculations the answer is clear, at large values of the external sources the classical solution is 
acceptable only at large values of the sources. In the symmetrical case of interactions with the small 
sources involved, the relative difference between the amplitudes is very large. The difference pretty 
fast decreases with the increasing of the value of the sources, but still the difference is not negligible 
even when q = 0.3 — 0.4, see Fig[H]-Fig[71 Important to underline, that one from the main reasons 
that provides the applicability of the classical solution is a presence of the third classical solution arose 
at some rapidity Yc ( see Section 2.2 an Eq. (j29p ) in the framework. The third solution provides the 
decrease of the classical amplitude at large values of the total rapidity and without this solution the 
precision of the model was even worse. The situation with the non-symmetrical case of interactions 
is better, see Fig[S]-Fig[21 Due the initially broken symmetry of the interaction, the fan diagrams 
are initially dominant and provide the better precision of the classical solution in comparison to the 
symmetrical case of interactions. 

Considering the second possible candidate on the role of "good amplitude approximation", eikon- 
alized function Eq. (j39|) . see Fig |12l we see that at small values of the external sources this amplitude 
is indeed better than the classical solution. In the region of large values of the sources both amplitude 
are equal more or less - the unitarization corrections are small for the large sources. We also could 
compare the effective Pomeron propagator with the asymptotic results obtained in the framework of 
s-channel model, see [35j . Looking in the Figllll and Fig[T6l we see, that the value of the full propa- 
gator, which considered as an amplitude in s-channel model, is close to one, especially for the second 
considered model with quaternary vertex included, which has direct relations with s-channel model, 
see again [35]. But, in general, for an arbitrary value of a source there is no propagator's unitarization 
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and the way to achieve the unitarization in this case is the eikonaUzation of the propagator only. 

Interesting problem, which also could be investigated in the framework of RFT-0, is a problem of 
the influence of the value of the triple Pomeron vertex on the behavior of the amplitude. Changing 
value of the p parameter from the 5 till 1 we change the value of the triple vertex A from 0.04 till 0.1. 
The resulting amplitude is depicted in the FigllOl From this plot we see, that already at relatively small 
values of rapidity this amplitude is decreasing till zero. This fast of the decreasing of the amplitude at 
large values of A denote a main difference between the models of RFT-0 with and without quaternary 
vertex. In RFT-0 without the quaternary vertex the amplitude, as it mentioned above, approaches 
zero whereas in RFT-0 with the quaternary vertex the amplitude approach constant with any value 
of triple vertex and correctly adjusted value of quaternary vertex, see Eq. (j44p . Eq. (j74p and Eq. (j76p . 
Of course, when the value of the triple vertex is small, two models do not so differ, see FigllSI This 
difference in the behavior of the amplitudes in two models may be also analyzed in the terms of 
s-channel model, see [53] for the details and explanations. 

As a possible example of the application of the model we calculated the values of the ratios of the 
single diffractive and differential single diffractive cross sections to the total cross section in the given 
framework, see EqlHT] and Eql88l The results of these calculations are present in the FigllSI - Fig|19l 
The plot Figll8l represents the ratio of single diffractive cross section (integrated over rapidity of the 
diffractive state) to the value of the total cross section. In two cases when the sources are equal 0.2 and 
0.7, we obtain, that for the fixed value of the diffractive state this ratio almost does not change with 
energy. More precisely, in both cases there is a tiny change in the ratio behavior at small rapidities 
4-5, and constant behavior at large rapidiies 6-8. Tracing the relative contribution of the diffractive 
state at fixed rapidity in both cases we see, that the contribution of the large mass diffractive state, 
Y2 = 4, approximately two times large than the contribution from small , "elastic", diffractive state 
I2 = 1, as it must be in reality. We see, that the simple RFT-0 model correctly reproduce main 
futures of the real QCD. It is also interesting to note, that the relative contribution of the large mass 
diffractive state is larger for the case of small external sources, there is R = 0.11 for source q = 0.2 
against R = 0.028 for source q = 0.7. Concerning the calculations of the second ratio, see Fig|19|. 
we obtain that these two cases are different. At small source, q = 0.2, the maximum of the ratio is at 
small diffractive state, which has almost constant behavior with total rapidity. At large vale of source, 
q = 0.7 , the situation is opposite, the maximim contribution comes from the large mass diffractive 
state. The contribution of the small diffractive state at g = 0.7 is zero in Figll9l at large values of Y, 
that means constant value of single diffraction cross section at small values of I2 and large values of Y . 
In general, considering Fig|18l - FigllDl we can conclude, that the description of the diffractive states 
in RET in the case when we include in calculations all possible re-scattering corrections is different 
from the "naive" diffraction models, where only part of corrections is included. Therefore, there is a 
hope, that using the same receipts of calculations in the QCD RET, we will be able correctly describe 
diffraction data and their energy dependence. 
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6 Conclusion 



In this paper we performed calculations of the different amplitudes in the framework of the RFT-0 
model. Comparing different approximation we analyze the applicability of the different approximation 
schemes and their dependence on the parameters of the model. The main conclusion of the calcula- 
tions is that we could trust to the classical approximation for the amplitude only when the value of 
the triple Pomeron vertex is small. In this region the difference between models with and without 
quaternary vertexes is negligible. Nevertheless, the important condition for this approximation in the 
case of interacting of symmetrical particles is an including of the third classical solution in the set of 
classical solutions, which arises at some critical value of rapidity. Without this solution the classical 
approximation is not good, even at not small values of the external sources. For the case of the small 
values of these vertexes the eikonalized Green's function amplitude is more precise solution then the 
classical one. Unfortunately, in real QCD this amplitude also could not be calculated precisely and 
therefore we could not trust to the classical solutions at all when the values of the external sources 
is not large enough. Does the proton is "large enough" in the case of real QCD is an open question, 
which could not be answered, unfortunately, in the given framework. 

When the value of the triple vertex is growing, the picture is changing drastically. In this case we 
could not trust anymore to the classical solutions of the models. Also, there is a large difference in 
the behavior of the amplitude in models with and without quaternary vertex, it's presence changes 
the asymptotic behavior of the amplitude. In real QCD it could means, that the different evolution 
equations must be applied in the different regions of the impact parameter space with the different 
values of the coupling constant. If we assume, that the influence of the non-perturbative effects is 
in the change of the value of the coupling constant only, then we must separate the contribution of 
the perturbative spots in impact parameter space, which evolution is described by the BK equation, 
from the non-perturbative regions where different evolution equations must be applied. In this case 
the overall amplitude is a sum of the amplitudes described by the different evolution equations in 
different regions, with the BK equation applicable in the regions of high partons density only. This 
picture will lead to the factorization of the non-perturbative effects from the perturbative ones, first 
of all, and also it could explain the possible applicability of the BK equation in case of proton-proton 
collisions. Namely, at high enough energies the overall contribution of "black" spots may be larger 
then the contribution of the "white" and "grey" ones and we could describe even inclusive data by the 
BK equation formalism. The large but constant non-perturbative contributions into the amplitude 
in this case may be accounted by the adjusting of the values of the external sources in interaction of 
interests. 

Concluding we underline, that RFT-0 model is a very interesting and important ground for an 
initial implementation of the different ideas which further could be applied in real QCD calculations 
as well. 
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